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ABSTRACT 

All prior work on measuring the spins of stellar-mass black holes via the X-ray continuum-fitting 
method has relied on the use of weakly-Comptonized spectra obtained in the thermal dominant state. 
Using a self-consistent Comptonization model, we show that one can analyze spectra that exhibit 
strong power-law components and obtain values of the inner disk radius, and hence spin, that are 
consistent with those obtained in the thermal dominant state. Specifically, we analyze many RXTE 
spectra of two black hole transients, H1743-322 and XTE J1550-564, and we demonstrate that the 
radius of the inner edge of the accretion disk remains constant to within a few percent as the strength 
of the Comptonizcd component increases by an order of magnitude, i.e., as the fraction of the thermal 
seed photons that are scattered approaches 25%. We conclude that the continuum-fitting method can 
be applied to a much wider body of data than previously thought possible, and to sources that have 
never been observed to enter the thermal dominant state (e.g., Cyg X— 1). 

Subject headings: accretion, accretion disks — black hole physics — stars: individual (H1743-322, 
XTE J 1550-564) — X-rays: binaries 



1. INTRODUCTION 

Black holes (BHs) are completely described by only 
three quantities: mass, charge, and spin. In astro- 
physical settings, any net charge will rapidly neutral- 
ize, with the result that a stellar-mass BH is specified 
by just its mass and spin. BH spin is commonly ex- 
pressed in terms of the dimensionless parameter a* = 
cJ/GM 2 with |a*| < 1, where M and J are respec- 
tively the BH mass and angular momentum, and c 
and G are the speed of light and Newton's constant. 
While mass measurements of stellar-mass BHs have been 
made for decades, the first spin measureme nts have been 
achieved only during the past three years (IShafee et al 
20061 : IMcClintock et all 2006; L iu et al.|[2008t IGou et a! 



200a Miller et ; al.ll2009l . and references therein). Mean- 



while, the spins of supermassive BHs have also been 
meas ured ([Brenneman fe Reynol ds 2006; iMiniutti et al.1 
2007). The only two methods presently available to 
measure BH spin are via modeling the thermal contin- 
uum spectrum of a BH accretion disk, as pioneered by 
I Zhang et al.l (|1997l ) , or by modeling the profile of a rela- 
tivistically broadened Fe K flu orescence line, as demon- 
strated bv lTanaka et all (|1995[ ). 

Spin is measured by estimating the inner radius of the 
accretion disk R lri . One identifies R m with the radius 
of the innermost stable circular orbit -RiscOj which is 
dictated by general relativity. i?isco /M is a monotonic 
function of a*, decreasing from 6 G?/c 2 to lG/c 2 as spin 
incre ases from a* = to a* = 1 ((Shapiro fe Teukolskvl 
1983). This relationship between a* and -ffisco is the 
foundation of both methods of measuring spin. 

In the continuum-fitting (CF) method, one determines 
-RiSCO by modeling the X-ray continuum spectrum, fo- 
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cusing on the thermal accretion-disk component. The 
observables are flux, temperature, distance D, inclina- 
tion i, and mass M. To obtain reliable values of spin, it 
is essential to have accurate estimates for M, i and D, 
which are typically derived from optical data. 

The CF method has been applied only to spectral data 
obtained in the thermal dominant (TP) state (or very re - 
cently to a near-TD intermediate state: IGou et al.ll2009h . 
The TD state is chiefly characterized by the dominance 
of the soft, thermal disk component of emission. (For a 
me asure of this dominance and a review of BH states, 
see lRemillard fe McClintoddl2006t ) The CF method has 
never been applied to the more Comptonized steep power 
law (SPL) state, which is characterized by the coexis- 
tence of a strong power-law component with photon in- 
dex r > 2.4 and a significant thermal component. Most 
models for the SPL state invoke Compton up-scattering 
of thermal seed photons by coronal electrons as the mech- 
anism that generates the power law. Herein, we employ 
a self-consistent Comptonized accretion-disk model that 
yields values of i?i n for SPL-state spectra that are con- 
sistent with those obtained for TD-state spectra. This 
result greatly increases the reach of the CF method, al- 
lowing one to obtain reliable measurements of spin for 
a much wider body of data than previously supposed, 
and for sources that do not enter the TD state (e.g., 
Cyg X-l). Moreover, the success of this model sup- 
ports the widely-held assumption that Comptonization is 
the mechanism which generates the observed high-energy 
power law component in SPL- and TD-state spectra. 

Our full model of a Comptonized accretion disk is a 
convolutio n of the relativ i stic thin accretion-disk m odel 
kerrbb2 (|Li et al.1 120051 : IMcClintock et all l2006t ) and 
SIMPL, an empirical model that convolves a Comptoniza- 
tion Green's function wit h an arbitrary seed photon spec- 
trum (jSteiner et al .1120091) . Both models are implemented 
in XSPEC (|Arnaudlll996h . SIMPL, with only two param- 
eters, ensures photon conservation and self-consistently 
generates the power-law component of the spectrum of 



an accreting BH using the accretion-disk component as 
input. 

We have chosen to apply our spectral model to the 
two bright transient X-ray sources H1743-322 (here- 
after H1743), which we feature, and XTE J1550- 
564 (hereafter J1550). Both transients are sources of 
large-scale relativistic jets an d high-frequency QPOs 
(jRemillard fc McClintockll2006l . and references therein). 
For a detailed comparison of the spectral and timing 
characteristics of these very similar transients durin g 
their principal outbursts, see iMcClintock et alj (|2009t ). 
Pres ently, the distanc e to J1550 is poorly constrained 
(see lOrosz et al.ll2002l ). and no useful distance estimate 
or dynamical information whatsoever is available for the 
black hole candidate H1743. Consequently, we cannot 
yet accurately estimate the spins of these black holes. In 
this work, we adopt fiducial values of M, i and D. Of 
course, i?i n (and a*) depend strongly on these fiducial 
values. However, as we show in H3.3I for any reasonable 
range of these input parameters, the dependence of R m 
on luminosity or on time during the outburst cycle is 
slight, which is an important conclusion of this work. 

We show that the very widely-used additive XSPEC 
models of Comptonization, namely the empir ical model 
POWE RLAW and the physical model COMPTT (|Titarchukl 
119941 : are inadequate for extracting measurements 
of spin from spectra with substantial power-law compo- 
nents. A self-consistent model such as SIMPL is required. 

2. OBSERVATIONS & ANALYSIS 

We apply the model described below to the full archive 
of spectral data for the 2003 outburst of H1743 (the 
most intense observed for this source) and for all five 
outburst cycles of J1550 obtained using the Rossi X-ray 
Timing E xplorer's (RX TE's) Proportional Counter Ar- 
ray fPCA: ISwankl ll999). We rely solely upon "standard 
2" spectra obtained using the PCU-2 module, RXTE's 
best-calibrated detector. All spectra have been binned 
into approximately half-day intervals, background sub- 
tracted, and have typical exposure times ~ 3000 s. For 
the first 5 weeks of PCA observations (through 2003 May 
1 UT) the detector was pointed 0.32° from H1743. We 
have corrected the fluxes to full collimator transmission 
assuming a triangular response with FWHM = 1°. We 
applied similar collimator corrections (« 0.1° — 0.3°) to 
three observations of J1550 performed on 1998 Septem- 
ber 7-8 and 1999 January 5 UT. 

A 1% systematic error has been included over all chan- 
nels to account for uncertainties in the response of the 
de tector (details on R XTE's calibration can be found 
in iJahoda et a l. 2006). As in our earlier work (e.g., 
IMcClintock et~aTl l200l . we have corrected for detector 
dead time while using contemporaneou s Crab observa- 
tions a nd the canonical Crab spectrum of lToor fc Sewardl 
(|1974h in order to calibrate the PCA effective area. The 
resultant pulse-height spectra are analyzed from 2.8 — 25 
keV using XSPEC vl2.5.0. 

In XSPEC, the model we employ is 
PHABS(simpl®kerrbb2), where pbabs is a widely-used 
model of low-energy photoabsorption. SIMPL redirects 
photons from the seed distribution, described here by 
the accretion-disk model kerrbb2, into a Compton 
power law. Like powerlaw, SIMPL has just two pa- 
rameters: (1) the fraction of seed photons fsc scattered 



into the power law, and (2) the photon power-law index 
r. SIMPL does not incorporate higher-order effects such 
as geometry-dependent scattering or reflection. The 
relativistic disk model kerrbb2 similarly has two fit 
parameters: (1) the spin parameter a*, which we express 
equivalently in terms of i?i n (33}, and (2) the mass ac- 
cretion rate M. From these two parameters we compute 
the Eddington-scaled disk luminosity, Lrj(a*, M)/ L^dd, 
where Lo is the luminosity of t he seed photons and 
L E dd ~ 1.3 x 1O 38 M/M erg s" 1 ([Shapiro fc Teukolskvl 
1983). The low -energy cutoff is parameterized in the 
PHABS component by the column density Ah, which 
we fix at 2 2xl0 22 cm~ 2 for H1743 and 8xl0 2 1 cm~ 2 
for J1550 (jMcClintock et all [20091 : iMiller et all [2001 . 
We include an additional model component to account 
for disk-reflection using the XSPEC model SMEDGE for 
J1550, which was not required for H1743. 

In our analyses described in § ^3.1I3.2\ we adopt the fol- 
lowing dynamical model parameters: For H1743, M = 
10 M Q , % = 60° and D = 9.5 kpc; and for J1550, 
M = 10 M Q , i — 70°, and D — 5 kpc. The values for 
H1743 are chosen arbitrarily to place the maximum out- 
burst disk luminosity at Ld / L-^dd ~ 0.7, and those for 
J 1 550 are round num bers based on the model described 
in lOrosz et alJ (|2002l ). In §3.31 we allow i and D to vary 
and consider six disparate dynamical models. 

For H1743 and J1550 we only select data over an order 
of magnitude in luminosity, between 0.05 < Ld/ L^dd < 
0.5 for the values of M, i, and D given above. This 
intermediate range in luminosity is chosen in order to 
eliminate both hard-state spectra that have little or no 
detectable thermal component and high-luminosity data 
for which the thin-disk approximation likely no longer 
applies. Further requiring goodness-of-fit (% 2 jv) < 2 and 
that the inner radius is well-determined (i?i n /Ai?j n > 
5, where Ai?j n is the ltr statistical uncertainty on i?i n ) 
leaves us with a total of 117 spectra for H1743 and 151 
spectra for J1550. 

We include for kerrbb2 the effects of limb darkening 
and returning radiation and set the torque at i?; n to zero 
(e.g. IMcClinto ck et al.l [20061 ). and for the dimensionless 
viscosity parameter we adopt a = 0.01. (Our results in 
the following section are modestly affected if one instead 
uses a = 0.1: R m is increased by ps 5% and becomes 
weakly dependent on luminosity, increasing by < 10% for 
an order of magnitude increase in Lr> . ) A color correction 
resulting from spectral hardening in the disk atmosphere 
is internally calcul ated for KERRBB2 using models KER- 
R BB and bbspec ([Davis fc Hubenvl l 2006h as described 
in IMcClintock et al.l ( 2006h . The upscattering-only im- 
plementation of SIMPL, known as SIMPL- 1, was used ex- 
clusively throughout unless otherwise noted. Larger val- 
ues of fsc are obtained using the double-sided scattering 
kernel SIMPL-2 (see Table 1), but i?j n and the other fit 
parameters are completely unaffected by the choice of 
kernel. 

3. RESULTS 

3.1. Final Selection of the Data via the Scattered 
Fraction 

The scattered fraction fsc sets the strength of the 
Compton power-law component relative to the disk. Fig- 
ure Q] shows for HI 743 the inner disk radius R[ n versus 
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TABLE 1 

Comparison of Bi n Across Spectral States 

Rin (in KERRBB2 used with 

















BH 


State 


N 


/sc a 


SIMPL-1 


POWERLAW 


compTT b 


H1743 


TD 


65 


0.012 


4.13 ± 0.05 


4.10 ± 0.06 


4.10 ± 0.07 




INT 


2 


0.062 


3.79 - 4.10 


3.48 - 4.02 


3.73 - 4.08 




SPL 


33 


0.135 


4.01 ± 0.06 


3.10 ± 0.24 


3.46 ± 0.27 


J1550 


TD 


100 


0.016 


5.20 ± 0.06 


5.05 ± 0.09 


5.14 ± 0.10 




INT 


18 


0.183 


5.16 ± 0.19 


4.37 ± 0.57 


4.93 ± 0.20 




SPL 


18 


0.123 


5.00 ± 0.15 


4.36 ± 0.21 


4.91 ± 0.26 



Fig. 1. — The inner disk radius i?i n versus the scattered fraction 
fsc for H1743. As indicated in the legend, the sy mbol types denote 
X-ray state ( see IRemillard fc Mc Clintock 2006). For / sc < 0.25, 
which is our adopted selection criterion, i?; n is generally constant; 
the median value for the TD-state data alone is indicated by the 
dashed line. However, for larger values of /sc> to the right of the 
vertical dashed line, the values of R ln diverge. Results are shown 
for all 117 spectra with x 2 l v < 2 and i?; n /Ai?,; n > 5 over the 
range Ljj = 5% — 50% I/Edd ( scc ©■ Error bars (lcr) on i?; n 
that are smaller than the plotting symbols have been omitted for 
clarity. Error bars on /sc are not shown; they are smaller than 
the plotting symbols except for extreme values of /sc (< 0.02 and 
> 0.6). 

/sc- For /sc < 0.25 the radius is quite stable and its 
value for the SPL data is very nearly the same as for the 
TD data. However, at large values of /sc the inner disk 
radius R m apparently recedes, indicating that either the 
model breaks down or a real change takes place in the 
dis k. One possible physic al explanation was proposed 
bv iDone fc Kubotal (|2006f ). who argue that in regimes 
of extremely-high Comptonization an inner disk corona 
can truncate the disk and increase i?j n by tens of percent, 
consistent with the high values shown in Figure [T] 

We have computed and compared plots of R ln versus 
/sc for four BH binaries (H1743, J1550, XTE J1655- 
40, and LMC X-3) and find that divergent behavior in 
their values of i?i n sets in for /sc > 0.2 — 0.3 (or fsc > 
0.25 — 0.4 for SIMPL-2). Based on a consideration of these 
results, we adopt /sc < 0.25 as a data-selection criterion 
in this work. The application of this criterion leaves a 
final data sample of 100 spectra for H1743 and 136 for 
J1550. 

3.2. Comparison with Other Comptonization Models 

Having applied our Comptonized accretion disk model 
phabs(simpl®kerrbb2) and obtained final data sam- 
ples for H1743 and J1550 ( §3.ip . we now analyze 
these selected data using alternative models for the 
Compton component. We employ (1) compTT, a 
widely-used model of Comptonization that describes 
the up-scatteri ng of blackbody- likc radiation by coro- 
nal electrons (iTitarchukl Il994f) . and (2) the empir- 
ical model POWERLAW. The full model formula- 
tions are respectively PHABS(kerrbb2+COMPTT) and 
PHABS(kerrbb2+powerlaw). We now use these mod- 
els to derive values of i?; n for both sources and compare 
these results to those obtained using our model. 

Figure [2] shows a side-by-side comparison of H1743 (left 
panels) and J1550 (right panels), where -R; n is now plot- 
ted versus Lu/L^dd (©■ The results in the upper pair 
of panels were obtained using our self-consistent Comp- 
tonization model SIMPL, while those in the lower pan- 
els were obtained using powerlaw. Plainly, for both 



Note. — The values and errors quoted for i?i n arc medians and 
median absolute deviations (MADs); we have chosen these quantities 
for their robustness. For Gaussian-distributed data, lcr ~ 1.5 MAD. 
i?i n is calculated using the fiducial Af , i, and D specified in ij2] 
a Calculated for SIMPL-1. For fits using SIMPL-2, /sc is »30% 
larger. k Geometry switch set to 1 (slab geometry) and redshift to 0. 
All other parameters are left free. 

sources SIMPL harmonizes the extreme discord between 
the SPL/intermediate (INT) data and the TD data that 
results from analyzing these data using powerlaw (Fig. 
[2f) & d). The reconciliation achieved using SIMPL (Fig. 
[2jz & c) indicates that the inner disk radii determined 
in the weakly-Comptonized TD state are very nearly 
the same as in the modcratcly-Comptonized INT and 
SPL states. Only data matching the selection criteria in 
§ ^213. II arc considered. 

Table [1] provides a summary of the results shown in 
Figure[2]and extends the comparison by including results 
for compTT. Qualitatively, the results for both sources 
are very similar; here we comment only on the results for 
H1743. Comparing SIMPL with POWERLAW, we see that 
for the former model R[ n is consistent between the TD 
and SPL states, 4.13 ± 0.05 and 4.01 ± 0.06, respectively 
(values and errors here are the median and median abso- 
lute deviation). On the other hand, POWERLAW delivers 
a radius for the SPL state that is « 24% smaller than 
for the TD state: 3.10 ± 0.24 versus 4.10 ± 0.06. While 
powerlaw fails dramatically to reconcile the TD- and 
SPL-state data, COMPTT provides only a modest im- 
provement, giving an w 16% smaller value of i?; n for the 
SPL state: 3.46 ± 0.27 versus 4.10 ± 0.07. The failure of 
COMPTT and POWERLAW to deliver a constant radius 
occurs because these additive models compete with the 
disk component for thermal flux and because they make 
no allowance for the flux which the disk contributes to 
the power law. 

3.3. Dependence on the Dynamical Model 

So far, our results are based on the specific and rather 
arbitrary dynamical model defined for each source in 
We now demonstrate that the quality of our results does 
not depend on the choice of a particular triplet of M, i, 
and D. For H1743 we analyze the data for six disparate 
dynamical models chosen as follows: We fix the mass 
at M = 10 M Q and vary the inclination from i = 30° to 
i = 80° in 10° increments, adjusting the distance in order 
to maintain the peak disk luminosity at Lu/Z/Edd ~ 0.7; 
this prescription leaves our selection criteria (§§2,3.1) 
largely unaffected. For this demonstration we restrict 
ourselves to a contiguous set of prist ine data that are 
free o f both edge and line features fsee lMcClintock et al.1 
[2009h . 
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Fig. 2.— The inner disk radius Ri n versus the Eddington-scaled disk luminosity Lo/LEdd f° r H1743 (left) and J1550 (right). Symbol 
types are defined in Fig. [T] For the upper pair of panels the Comp tonization model employed is SIMPL and for the lower panels it is 
POWERLAW. The data sample considered here is that described in EI3.1I For J1550 note in panel c the many INT-state data that are brought 
into agreement with the SPL- and TD-state data when applying SIMPL. Error bars are omitted when smaller than the symbols. 

Figure [3Jj shows a portion of the 2003 outburst light 
curve of H1743. Figure[3}? shows corresponding values of 
Rin versus time for the six models described above. We 
draw the following key conclusions from Figured) (1) R ln 
is constant for each model to within s=s 2% as the source 
passes from the SPL state to the TD state, and as the 
source flux decays by a factor of w 6. We furthermore 
note that R[ n is stable during the two strong SPL-state 
flares that occur on days 75.6 and 79.5. (2) The character 
of the small systematic variations that occur in i?j n dur- 
ing this entire 4-month period are essentially the same 
for all six models. For completeness, we recomputed all 
the results shown in Figure \3}> using first M = 5 M and 
then M = 15 M Q . Apart from offsetting the value of 
i?in, the character of these results is the same, including 
the level of scatter, as for the case of M = 10 M Q . We 
conclude that, apart from setting the median value of 
Ri n , the choice of model has no significant effect on the 
results presented in Figure [3] 

Likewise, for J1550 we analyzed a ~ 130-day stretch 
of data obtaine d during the 1998 o utburst cycle (MJD 
51110 - 51242: ISobczak et al.ll2000fh We assumed fidu- 
cial values of M and i and explored a wide range of dis- 
tances from D = 3 — 8 kpc. We obtained results very 
similar to those presented for H1743 (Figure \5f>) , consis- 
tent with an internal scatter of w 2%. 
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Fig. 3. — (a) A contiguous 126-day portion of the 225-day RXTE 
PCA light curve of H17 43, which is shown in full in Figure 3a of 
IMcClintock et all 120091) . The 2-20 keV unabsorbed fluxes were 
obtained by modeling the PCA spectral data. Time zero is the date 
of discovery of H1743 during its 2003 outburst, which occurred on 
2003 March 21 (MJD 52719). (6) R in versus time for the six models 
described in the text, shown as alternating black/blue tracks for 
clarity. The median absolute deviations for the extreme models 
with t = 30° and i = 80° are 2.2% and 1.8%, respectively. Fluxes 
for 83 spectra are plotted in panel a and 79 values of R ln are plotted 
in panel b (except for i = 40° with 78); i.e., four (five for i = 40°) 
spectra failed to meet our selection criteria. Error bars are omitted 
where they are smaller than the symbols. 



4. DISCUSSION 



iKubota et all (|200lD and lKubota fc Makishimal (|2004T ) 
present the first self-consistent treatment of disk- 
dominated accretion at high luminosity in black hole bi- 
naries. They showed for GRO J1655-40 and J1550 that 
what previously had appeared to be anomalous behavior 
was a natural result of strong inverse-Compton scatter- 
ing. In particular, they demonstrated that the inner disk 
radius was stable when the flux attributed to the power 
law was properly associated wit h the disk. Their res ults 
have been confirmed recently bv lSteiner et al.1 (|2009h us- 
ing SIMPL (§§1,2). In this paper, we provide additional 
support for the work of Kubota et al., while supplying 
in this context the first relativistic analysis of the accre- 
tion disk component. Both the earlier work by Kubota 
et al. and this Letter demonstrate that, when modeling 



Comptonization, a self-consistent treatment is necessary 
in order to explain BH behavior across spectral states. 

In all of our earlier work measuring the spins of BHs 
using kerrbb2, we have selected data with L£)/LEdd < 
0.3, which corresponds to t he thin-disk limit (H/R < 
0.1; IMcClintock et all [20061 ). In the present work, the 
luminosity of J1550 is very uncertain and that of H1743 is 
unconstrained. For this reason, we present a broad range 
of luminosities, which likely exceeds the thin-disk limit. 
In work aimed at determining BH spin, when reliable 
distance estimates and dynamical data are available, one 
should apply the aforementioned luminosity restriction. 

In conclusion, we have analyzed a selected sample of 
~ 100 spectra for each of two bright transient sources 
using the self-consistent Comptonization model SIMPL 



4 



convolved with a rclativistic accretion disk model. We 
have thereby shown that the derived inner disk radii - or, 
equivalently, the derived spins of these BHs - remain sta- 
ble to a few percent whether the source is in the TD state 
or the more strongly-Comptonized SPL state. We have 
further shown that this stability holds for fsc ^ 0.25 and 
for a wide range of input model parameters. We conclude 
that the continuum-fitting method of estimating BH spin 
can be applied to far more X-ray spectral data and more 
sources than previously thought possible. 
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